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Q_i! We present a simple model with string absorption and percolation to describe 

o ■ the J/ip suppression in heavy ion collisions. The NA50 data are fairly well 

■ explained by the model. 
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It was recently argued that central distributions, in transverse energy or multiplic- 
ity n, are dominated by fluctuations in the number v of elementary collisions, in the case 
of nucleus-nucleus interactions and even hadron-hadron interactions at very high energy 
0- 

In that framework, if the total number of events with v elementary collisions is N(v) 
and Nc(v) is the number of events with the rare event C occurring, then |fj| 

N c {y) = a c vN(v), (1) 

where ac is the probability of event C occurring in an elementary collision. An event is 
rare if, in good approximation, it only occurs once. 

Relation ([]]) is correct as far as event C is free from absorption 0, or from quark- 
gluon plasma formation 0. This is directly seen in the observed linear relation between 
E T distributions associated to Drell-Yan productions, Nd.y,(Et), and to minimum bias 
events, N(E T ) g |. 
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However, in the case of J/0 production, the formed J/0 may be destroyed by subse- 
quent interactions, as it moves through the interacting medium ||, or it may happen that 
it is not really formed, as a result of quark-gluon plasma Debye screening preventing the 
cc binding. 

In the standard multicollision model — dual parton model, 0, for instance — in every 
elementary collision two coloured strings are originated. These strings constitute the 
medium that may absorb a created J/0 (see also ||). They may also fuse (percolate) and 
form the quark-gluon plasma, thus inhibiting J/-0 creation. 

We shall consider that when v collisions occur they occupy an interaction volume Viy) 
characterised by a transverse area A{y) and a mean longitudinal length L(v), 

V{v)=A{ V )L{u). (2) 

As nuclear matter density is uniform, we expect the string density p s , 

2v 

Ps = A{u)L{uY (3) 

where 2v is the number of formed strings, to be as well uniform. Note that p s in eq. (|3|) 
is, in general, independent of the interacting nuclei, but may be energy dependent. 

In order to take care of J/-0 absorption we introduce the survival probability P s {v) in 
the conventional form 

P s (v) = exp(-L(v)p a o), (4) 

a being the J/-0-string (qq) absorption cross-section. 

Next, we include the probability of quark-gluon plasma formation in a very simple 
two-dimensional percolation model |9|, |j~Q| . In a v elementary collision event there is a 
probability P nv {y) of percolation not occurring and a probability 1 — P nv {y) of percolation 
occurring. As percolation means here quark-gluon plasma formation, we make the strong 
assumption that the J/0 is formed only in events in which there is no percolation. This 
means that in eq. ([!]) we shall multiply Nc{v) by P 7 



np 



v): 



N c (u)-^N c (u)P np (v). (5) 

This probability P np (^) is in fact a function of the scaling variable rj, the dimensionless 
density of strings in the transverse plane, 

2 2u 

^-r s — , (6) 

where nrl is the transverse area of a string. Because of eq. (|3|), we can also write 

r] = nr 2 sPs L(v). (7) 



Recently, in [11 , the probability P np was studied for S-U and Pb-Pb collisions at 



a/s = 19.4 AGeV. It can be written in the form 
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Figure 1: Percolation probability as a function of the transverse density rj for S-U (open squares) and 
Pb-Pb (filled circles) collisions geometry as obtained from Monte Carlo simulations |ll]]. The curves 



show fits with P p {ri) = 1 — P np (rj), see eq. (]8|) ; we obtain a 
a = 0.04 ± 0.01 for Pb-Pb (solid line); 



= 0.07 ± 0.01 for S-U (dashed line) and 
r] c = 1.15 ± 0.02 



where r\ c is the critical density, r\ c m 1.15 for a uniform distribution of strings in the 
transverse plane, and a is a parameter which depends only on the geometry of the colliding 
nuclei (see Figure p. 

We can now write our final formula for production, normalised to Drell-Yan 
production, R = (J/i/>)/(D.Y.), as 

R{u) = KP s (v)P np (is) (9) 

where K is a normalisation constant, which should be close to 50 from pp collisions |1J 
This corresponds to the small v, small Et limit, with P s (v) — > 1, Pnp{v) — * 1- Validity of 
eq. (H) requires J/?/> and Drell-Yan production being treated as rare events. 

Now we note that, from eq. (f|) and eq. (0), we obtain that R is in fact a function of 
the mean length L only, 



R = K exp(—Lp s a) 



-i -i 



(10) 



For small values of L, away from percolation, the ratio R in eq. (|10|) is universal, i.e., does 
not depend on the colliding nuclei, but it depends, in our approach, on the cm. energy, 
through p s . 

Results were presented by the NA50 Collaboration || on S-U and Pb-Pb collisions as 



functions of the average interaction distance L computed in a geometrical model ||T3j, |T4 
In Figure ||a we compare eq. (|10|) with the NA50 data. The string radius was fixed at the 
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Figure 2: The ratio of 3/ipto Drell-Yan events as predicted by eq. ( 

as published by NA50 Collaboration jB 



compared to experimental data 



value r s = 0.2 fm, |TD|, which places r\ for S-U somewhat below the percolation threshold, 
and the parameter a at the Pb-Pb value, 0.04. In Figure ||a we have included the 1998 
data on the ratio R as a function of E Tl and transformed to L using the NA50 E T vs L 
plot |13| in the approximate form (valid for Pb-Pb at 19.4 AGeV), 



L « 2.07^^, 



ill) 



with Et in units of GeV. 

The curve shown in Figure ||a corresponds to the following values of the parameters 
in eq. (|i~0|): p s = 0.9 fm -3 , a = 1.5 mb and K = 58. The string density p s is, as expected 



15fl , much larger than nuclear density (p ~ 0.17 fm~ a ) and the absorption cross section 
agrees with what is measured in J/?/;-hadron collisions at intermediate energies. 

In Figure Qb we present the ratio R as a function of Et in comparison with data, by 
making use of eq. (ITT]) . The agreement is mostly qualitative. 

The essential point is that at small E T the behaviour is typical of absorption — with 
positive curvature — at higher values of E? the approach to the percolation transition 
produces a change in the sign of the curvature. 

Finally, we present in Figure |3| our predictions for RHIC and LHC energies, compared 
with the curve for SPS. Figure ||a displays eq. ([10| ) with the expected change in p s deter- 
mined by the change in the number of strings given in ||10||, and with the other parameters 



(a, r s ) held constant, In Figure [|b we transformed the curve from the L variable to Et 
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Figure 3: The ratio of 3/tp to Drell-Yan events as predicted by eq. (0), predicted at RHIC energy (dashed 
lines) and at LHC energy (dotted lines), compared with SPS energy (solid lines). 



with the conservative assumption that the proportionality between Er and the number v 
of elementary collision is independent of yfs. This implies that 



(12) 



and we use eq. (|TTj) for y/so = 19.4 GeV. We see that, as ^fs grows, both the absorptive 
part gets larger and the change of curvature denoting percolation sets in at a much smaller 
value of E T . 
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